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Abstract 



The Seyfert 1 galaxy, NGC3516 has been the subject of many absorption line studies at 
both ultraviolet and X-ray wavelengths. In the UV, strong, broad, variable associated metal 
line absorption with velocity width ~ 2000 km s _1 is thought to originate in gas with Nh 
>10 19 cm -2 lying between 0.01 and 9 pc from the central active nucleus. The Ginga X-ray 
data are consistent with several possibilities: a warm absorber and a cold absorber combined 
either with partial covering or an unusually strong reflection spectrum. We present ROSAT 
observations of NGC3516 which show a strong detection of a warm absorber dominated by 
a blend of OVII/OVIII edges at ~ 0.8 keV with N H ~7x 10 21 cm" 2 and U: 8-12. We 
argue that NGC3516 contains an outflowing 'XUV absorber showing the presence of X-ray 
absorption edges, that are consistent with the presence of broad absorption lines in the old 
IUE spectra and their disappearance in the new UV observations. Our dynamical model 
suggests that the OVII absorption edge will continue to weaken compared to the OVIII 
edge, an easily testable prediction with future missions like AXAF. Eventually the source 
would be transparent to the X-rays unless a new absorption system is produced. 
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1. Introduction 

Given the recent success of efforts to unify associated ultraviolet (UV) metal line 
absorbers with ionized ("warm") absorbers seen in the X-ray (Mathur et al. 1994, 1995, 
Mathur 1994), NGC3516 is a prime candidate for such an investigation. NGC 3516 is a 
low-redshift (z=0.0089) Seyfert 1 galaxy with a variable, associated UV absorption line 
system which has been extensively monitored with IUE (Voit et al 1987, Walter et al 
1990, Kolman et al. 1993 (KEA); Koratkar et al. 1996). Based on IUE spectra, NGC 3516 
contains the strongest (largest equivalent width) UV absorption system known in a Seyfert 
1 galaxy. This system contains at least two distinct components: a broad (FWHM~2000 
km s^ 1 ) variable component, most likely to be associated with the X-ray absorber, and a 
narrow (~ 500 km s" 1 ) non-variable component. The equivalent-width of the CIV A1549 
broad+narrow absorption line varies by a factor of 4-5 and anti-correlates with the UV 
continuum level on short timescales of days. Walter et al. (1990) have also found variations 
on longer timescales (~ years) possibly due to the absorbing cloud passing through the line 
of sight. Recent observations have found that the broad variable absorption lines of CIV 
and NV have disappeared (Koratkar et al. 1996, Kriss et al. 1996b). 

The X-ray absorber is not so well-known or heavily studied. The highest quality 
data to date are from 1989 Ginga observations reported by KEA, consistent with several 
possibilities: a warm absorber, a cold absorber with partial covering, or a cold absorber with 
an unusually strong reflection spectrum. If the absorption is due to the cold matter, the 
column density is large N# ~ 10 22 cm -2 . These authors concluded that a common origin 
for the X-ray and UV absorbing material is not possible due to the apparently different 
ionization states and column densities of the UV and X-ray absorbers. Ginga data were 
later re-analyzed by Nandra & Pounds (1994), showing that a warm absorber describes the 
data well. New ASCA observations of NGC3516 are presented in a recent preprint by Kriss 
et al. (1996a), which show a presence of an ionized absorber similar to the one discussed in 
this paper. 

Here we present a long (~ 1 day) high signal-to-noise ROSAT (Trumper 1983) Position 
Sensitive Proportional Counter (PSPC, Pfefferman et al. 1987) spectrum of NGC3516 
obtained in 1992 which, with its lower energy range, is sensitive to absorption edges due to 
ionized oxygen. The high S /N allows the strong detection of both OVII and OVIII edges 
independently, in spite of the limited spectral resolution of the PSPC. These edges show 
the presence of an ionized absorber in X-rays. We discuss the constraints these data place 
on the physical conditions in the absorber and revisit the question of a common origin for 
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the X-ray and UV absorbing material. 



2. ROSAT Observations and Data Analysis 

NGC 3516 was observed with the ROSAT PSPC on October 1, 1992 for a total 
livetime of 13,081 sec. over a real time span of 83,455 sec. We retrieved the data from 
the HEASARC0 database and analyzed it using the PROS0 package in IRAF. The source 
counts were extracted from within a 3' radius circle centered on the source centroid. The 
small background was estimated from an annulus centered on the source and inner and 
outer radii of 5' and 7' respectively. The total net counts were 59,225±246 yielding a count 
rate of 4.53±0.02 s _1 (Table 1) (c.f. 1420 ± 340 counts for 3C351 in which absorption edges 
were detected, Fiore et al. 1993). 

Since NGC 3516 is known to be a variable X-ray source (Halpern 1982), we looked for 
X-ray flux variations over the total span of the observation. A total of eight good time 
intervals (GTI) were found. The net counts in each GTI were extracted in the same way 
as discussed above. The live time, count rate and the total counts in each GTI are given 
in Table 1. The change in count rate shows the flux variability which is clearly seen in the 
light curve in Figure 1. The variations are significant at the 99.99% level (x 2 test, K-S test 
using PROS task "vartst"). Since variability of ~40% is clearly seen, the source spectrum 
for each GTI was extracted separately, and analyzed using XSPEC0. As can be seen in 
Table 1, the individual GTIs contain sufficient counts to fit complex spectra. 



2.1. Spectral Analysis 

The results of the spectral analysis of all the datasets are presented in Table 2. All 
the spectral fits were made to the 3-34 PHA channels as extracted by the standard PROS 
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analysis. Channels 1 and 2 were ignored since they are inadequately calibrated. The 
response matrix released in January 1993 was employed. The errors represent the 90% 
confidence interval. 

Our first fit used a simple power law with absorption by cold material at zero redshift 
with solar abundances, and cross sections by Morrison and McCammon (1983). Both the 
power-law slope and absorption column density were free to vary. In all GTIs, the column 
density of the cold absorber (N#) was consistent (within 90% confidence) with the Galactic 
value of 3 x 10 20 atoms cm~ 2 (Heiles 1975) (e.g. = 3.0 ± 0.1 x 10 20 cm" 2 in GTI 4; and 
Nh = 2.4 ± 0.5 x 10 20 cm~ 2 in GTI 7). The single power-law (and all subsequent) fits were 
redone with N# fixed at the Galactic value and the results are presented in Table 2. 

The single power-law fits are unacceptably bad, with xl ~4 and significant negative 
residuals between 0.7 and 1 keV (see Figure 2). This is a clear signature of K-shell 
absorption edges due to OVII and/or OVIII. We then fitted the spectra with a power-law 
(PL), cold Galactic absorption, and an edge at the energy of a redshifted OVII K-edge 
(E(rest frame)=0.74 keV, z=0.0089). The improvement in the fit was dramatic, as 
demonstrated by the results of an F-test (Table 2), although small residuals remained 
near the edge in several GTIs. For consistency we then fitted the spectrum of all GTIs 
with an additional edge at the energy of the redshifted OVIII K-edge (E(rest frame)=0.87 
keV). Once again the improvement in the fit was significant and the spectrum was well fit 
(xl — 1-1) Table 2, Figure 3) except for GTI 4 (xl = 1.6, see below). Thus, even though 
the energy resolution of the ROSAT PSPC is insufficient to resolve the OVII/OVIII edges, 
the absorption in NGC 3516 is so strong (r ~ 1) that the data require two separate edges 
in half of the GTIs (> 99% confidence, labeled " c " in Table 2). These are not simply the 
four GTIs with highest S/N. 

As discussed above, the data were consistent with N^=N^ (Galactic). To confirm that 
this is the case even with the 'power-law plus two edges' model, we fit the spectra with a PL, 
Galactic N# (fixed), two edges and an additional cold at the source. The fitted value 
of the additional column density was always much smaller than Galactic (e.g. 1.2xl0 18 
atoms cur 2 for GTI 1) and consistent with zero. Excess is clearly not required. GTI 4 
was the only exception for which an excess absorption N# = 4.9 ± 2.6 x 10 19 improved the 
fit significantly (xl = 1.28, improvement > 99% significant, F-test). Parameters for this 
best fit model are given in the last line for GTI 4 in Table 2. Since GTI 4 has the maximum 
S/N, it is possible that a small additional column is present (Nfl-<3.7 x 10 19 cm~ 2 ) but 
undetected in the remaining GTIs. 
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Since the flux is found to be variable but the absorption and spectral shapes are not 
(see section 3.1), we also fitted all the GTIs together with a 'power-law plus two edges' 
model allowing their normalizations to be free (Table 2). Excess absorption was not 
required, though the best value was N# = 3.3^;!>;3 x 10 19 cm -2 , similar to the results above. 



3. The X-ray warm absorber 

Knowing the total column density and the opacity of the oxygen edges, we can 
constrain the ionization state of the warm absorber (Mathur, Elvis & Wilkes, 1995). The 
present case of NGC3516 is a little tricky, however, because N# is not known. Instead we 
need to estimate the N# and the ionization parameter (U) by fitting the spectrum with a 
warm absorber model. Alternatively we can derive these parameters in a more elegant way, 
exploiting the fact that we observe both OVII and OVIII edges in NGC3516. We apply both 
these methods below. 

Let tqvh and tqviii be the opacities due to K-edges of OVII and OVIII respectively. 
The absorption cross-sections of OVII and OVIII are 0.28 xl(T 18 cm" 2 and 0.098 xl0~ 18 
cm -2 respectively (CLOUDY; Ferland 1991). Using these we can derive the column 
densities of the two ions, Nov// and Now// and the ratio ^ovii/^oviii = iovii/iovin 
where iovn and iovin represent the fraction of oxygen in these two stages of ionization. 
Over the span of ROSAT observations (GTI 1 through 8), iovn/iovin ranges from ~ 0.16 
to 5.8. We determined the dependence of iovn and iovin on the ionization parameter, 
U, for a photoionized gas cloud using CLOUDY and assuming solar abundances and a 
standard AGN continuum (Ferland 1991, version 80.06). The density was assumed to be 
10 5 cm~ 3 . The results are displayed in Figure 4 and demonstrate that the observed range 
in iovn/iovin occupies a small area in the parameter space of fractional ionization as 
a function of ionization parameter. Using Figure 4, we can not only infer the ionization 
parameter (U ranges from ~ 3.2 at GTI 6 to ~ 18 at GTI 3), but also iovn and iovin- 
These lie in the range 0.06< fovn <0.8, and 0.1< fovm <0.5. Using the abundance of 
oxygen relative to hydrogen (8.51 xlO -4 ; Grevesse & Andres 1989) and iovn or iovin-, the 
total column density is readily calculated to be N# = 0.4 — 2.2 x 10 22 cm -2 . 

As discussed above, the ionization parameter of the X-ray absorber is determined using 
Figure 4. It should be noted however that the exact relation of ionization fractions and U 
depends upon the shape of the input continuum (Fiore et al. 1993, Mathur et al. 1994). We 



- 6- 



have used a standard AGN continuum (CLOUDY, Ferland 1991) in our analysis. With a 
different continuum shape, the inferred range of U would change but, qualitatively, Figure 
4 would remain the same. For example, an X-ray slope of 1.0 (as seen in the ROSAT data) 
rather than the 0.7 assumed in the standard continuum leads to a larger value of U: 8.5 < 
U t 23. On the other hand if the X-ray slope is much flatter (a ~ 0.3) as seen by Ginga, 
U would be correspondingly lower. Independently of the shape of the input continuum, 
the X-ray data for the absorber in NGC3516 requires iovii/iovm ~ 1) corresponding to 
a highly ionized absorber. A similar result was obtained by Kriss et al. (1996a) who used 
a variety of incident continua and found that their results were independent of the exact 
shape of the continuum apart from the deduced ionization parameter U. The only derived 
quantity that depends upon the exact value of U (and so on the shape of the continuum) 
is the distance of the absorber from the central continuum source (R 2 = Q/AnUnc). Since 
there is a four orders of magnitude uncertainty in density, n (section 3.1), the uncertainty 
in R is very large, far exceeding that due to the uncertainty in U. We will thus assume the 
standard continuum in the rest of the paper. Kriss et al. (1996a) use a continuum with 
significantly flatter a ox and correspondingly lower value of U=1.66 in their analysis of 1995 
ASCA data, so we cannot directly compare our U value with Kriss et al. (1996a). 

Unlike the other objects showing an X-ray warm absorber (e.g. 3C351: Fiore 
et al. 1993, NGC3783: Turner et al. 1994, NGC5548: Mathur et al. 1995), the X-ray 
continuum in NGC3516 is transparent at low energies. To investigate this absence of a 
low energy turnover in the ROSAT spectrum of NGC3516, we generated the transmitted 
spectrum for the best fit U and N# (using CLOUDY). Figure 5 shows the input spectrum 
(solid line; "standard" AGN spectrum) and the transmitted spectrum (dotted line) over the 
ROSAT and Ginga energy bands. It is clearly seen that the spectrum is transparent at low 
energies and the OVII/OVIII absorption edges are prominent in the ROSAT energy range. 

As a consistency check we also fitted the spectra with a warm absorber model. The 
models were generated using CLOUDY and then incorporated into XSPEC as 'table 
models'. The input parameters to CLOUDY were same as used above. The results of the 
fits are given in Table 3 and are consistent with the results discussed above. When N# is 
fixed to the Galactic value, the fits are acceptable for all the GTIs except GTI 4, similar to 
the results in section 2.1 (first line for each GTI in Table 3). When N# was allowed to be 
free, the fits improved and were good in all the GTIs (second line in Table 3). The best fit 
N# was slightly larger than but consistent with the Galactic except in GTI 4, again similar 
to the results above. In these models the range of the ionization parameter is 3 < U < 28 
and the column density of the warm absorber is N#=0.4-3.2xl0 22 cm -2 , consistent with 
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the results of the empirical fits. 



3.1. Variability 

As mentioned previously, the flux level of NGC3516 is variable as a function of time 
(Figure 1). However, the variability of opacities in OVII and OVIII absorption edges is not 
significant (Figure 1, <2<r variations). So all the GTIs were fitted together, leaving relative 
normalization free, resulting in a tighter constraint on U (6^U^13, 90% confidence). If the 
entire data set is fitted (i.e. without dividing into separate GTIs), then the range in U is 
even smaller (7.9<U< 12.6, 90% confidence) and N^=0.7±0.1 x 10 22 cm" 2 (Table 3). We 
use these values of U and N# in the rest of the paper. 

The ionization parameter is expected to increase with increasing flux on the 
photoionization time scale. The OVII photoionization time scale is t ph = 2 x 10 5 ng 1 seconds 
(see e.g. Reynolds et al. 1995). The lack of variability in OVII the opacity implies that 
tphZW 1 seconds. This puts an upper limit on the density of the absorber, n<2 x 10 7 cm -3 . 

In the UV, variations of absorption line strengths on timescales of days have been 
discussed for the CIV absorber and used to derive a density of n> 10 5 cur 3 in the absorbing 
region (Voit et al. 1987). This calculation, however, assumes that CIV is the dominant 
ionization state. In the present case, a correction factor of nciv/ n cv needs to be applied 
to the recombination time scale. The resulting lower limit on the density of the absorber is 
then smaller, n>10 3 cm~ 3 . 



4. Comparison with the Ginga results 

The power-law slope in the ROSAT PSPC data is much steeper (a ~ 1) than that seen 
by Ginga (0.25 < a < 0.43, KEA). Since the PSPC range is dominated by the absorber, 
the Ginga slope is more likely to indicate the true power-law slope. The ROSAT PSPC 
data show no significant excess cold absorption, but require a warm absorber with N H 
~7x 10 21 cm -2 and 8<U<13. It can be seen from Figure 5 that this model does not 
predict the presence of an Fe edge (r pred < 0.025). This is inconsistent with the Ginga data 
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which showed a "cold" absorber (N H ~ 10 22 cm -2 ) and an Fe edge. As can be seen in 
Figure 5, the ionized absorber causes a low energy turn-over within the Ginga band which 
would mimic "cold absorption" in the Ginga data. Indeed, re-analysis of the Ginga data 
showed that the low energy turn-over is most likely due to the warm component (Nandra 
& Pounds 1994). However, the origin of the Fe-K edge would have to be different from the 
XUV absorber discussed here. It may be, e.g., from the torus (Krolik and Kriss 1995) which 
possibly grazes our line of sight. Alternatively, the Fe edge and the large column may be 
associated with an absorption system that has moved out of the line of sight. We note in 
passing that there is no obvious Fe edge in recent ASCA observations (Kriss et at. 1996a). 



5. Comparison of X-ray and UV absorbers 

In many AGN, the X-ray and UV absorbers have been found to be one and the same 
(the 'XUV Absorbers': Mathur et al. 1994, 1995). In these cases apparent inconsistencies 
between the conditions in the UV and X-ray were resolved by realizing that the ions 
observed in the UV were not the dominant ions in the gas. NGC 3516 also shows both X-ray 
and UV absorbers. In a previous study based on the simultaneous (October 1989) Ginga 
and IUE observations, KEA concluded that, if an X-ray ionized absorber is present, it is 
unlikely to be due to the same gas as the UV absorber. The combination of the strength 
of the UV CIV absorption and the very high ionization state of Fe (>12, consistent with 
> 5x ionized C) indicated by the ionized absorber model fit to the Ginga data implied to 
the authors that a consistent solution was unlikely, although no detailed calculations were 
presented. Since the ionized absorber observed with ROSAT is clearly different from the 
Ginga absorber, the question of an XUV absorber in NGC3516 should be revisited. 

Here we investigate quantitatively whether the ROSAT ionized absorber is consistent 
with the broad UV absorber with high ionization lines. We use the method described 
in Mathur et al. (1995). The ROSAT observations were made in 1992. There were no 
simultaneous UV observations. The 1993 IUE observations show that the broad UV high 
ionization absorption lines had disappeared (Koratkar et al. 1996). Figure 4 shows the 
model with fractional ionization of CIV as a function of U. If the ionization parameter of 
the X-ray absorber were toward the high end of the observed range (U=12.8), the ionization 
fractions of CIV and NV would be very small (log f 4). Assuming solar abundances of 
carbon and nitrogen (3.63 xlCT 4 and 1.12 x 1CT 4 respectively, Grevesse & Andres 1989) 
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and the total column density of 7 x 10 21 cm~ 2 as derived from the X-ray data gives N(CIV) 
< 10 14 cm~ 2 and N(NV) < 8 x 10 13 cm -2 . The absorption lines would not be detected. 
In this case there would be an X-ray ionized absorber, but no broad high ionization UV 
lines of CIV or NV, consistent with the observations. However, the X-ray absorber must 
produce a detectable OVI A1034 line with — 2.1 < log iovi ~ — 2.7 (Figure 4). Since there 
were no contemporaneous far-UV observations with ROSAT, this cannot be directly tested. 
Analysis of absorption lines embedded in the emission line profiles of quasars is extremely 
difficult since the shape of the emission lines is unknown. However, visual inspection of 
the 1995 HUT data for NGC3516 (Kriss etal 1996b) shows a broad absorption feature 
blue-wards of the OVI emission line peak and at roughly the expected redshift of the broad 
UV absorber. 

It is also of interest to point out the position of the broad CIV and NV absorption lines 
on the f-U curve when they were strong. The equivalent width of the variable broad CIV 
absorption line in the IUE observations ranged from ~3A to ~10A (KEA) leading to a 
lower limit on the CIV column density Nciv > 2.2 x 10 15 cm -2 . Similarly Voit et al. (1987) 
report Nqiv ~10 15 cm -2 . Assuming this lower limit gives a lower limit on the ionization 
fraction of CIV: W~4 x 1(T 4 (log f C /y~-3.4). The NV EW given by KEA varies in the 
range 1.1 - 7.1 A. This corresponds to N(NV) > 3.5 x 10 14 cm 2 (using oscillator strength 
f=0.235, Allen 1973). The lower limit on the ionization fraction of NV is estimated to be 
f7vy~5 x 1CT 4 (log fAry>-3.3). The arrows in Figure 4 indicate these lower limits on i C iv 
and i NV observed by IUE. These lie in a range of U smaller than that corresponding to 
the ROSAT observations. This suggests that if U increases as a function of time, both the 
former presence and current absence of the broad UV lines would be consistent with the 
observed X-ray warm absorber. 

The IUE data have occasionally shown the presence of a noisy, possibly broad SilV 
absorption feature (Voit, Shull & Begelman 1987). The conditions of the warm absorber 
derived here do not predict detectable SilV absorption (log fsuv~— 30). However, SilV 
absorption is also present, with no significant change in strength, in the 1993 IUE spectra 
in which no broad, high ionization lines are reported (Koratkar et al. 1996). Thus it 
appears that the feature is a part of the narrow UV absorption system rather than the XUV 
absorber. Note that the two other objects with well-studied XUV absorbers (e.g.3C351, 
NGC5548) do not have any such additional narrow absorption systems. It is possible that 
the second, low column density absorber observed in the ROSAT data is related to this 
narrow UV system. 



6. Discussion 



We argue that NGC3516 contains an XUV absorber. It has high column density 
(N# = 7 ± 1 x 10 21 cur 2 ) derived from the X-ray observations. The UV observations 
imply that the absorber is outflowing with a velocity of ~ 500 km s _1 (given its blueshift) 
and the line width implies an internal velocity dispersion of < 2000 km s _1 (see sec. 1). 
Assuming the upper limit on its ionization parameter (U<13), a density of 2 x 10 5 cm~ 3 , 
and N# = 7 x 10 21 cm~ 2 , the distance of the absorber from the central continuum source is 
< 5 x 10 18 cm and its thickness is about 3.5 x 10 15 cm. At this high value of U no detectable 
broad UV lines are predicted (sec. 3, fig. 4), consistent with the disappearance of the broad 
UV absorption lines ( Koratkar et al. 1996, Kriss et al. 1996b). 

Kriss et al. (1996b) have presented higher resolution ASCA data of the ionized absorber 
in NGC3516. Based on the non-detection of X-ray absorption lines in the ASCA data 
they have argued that the velocity dispersion parameter 'b' of the absorber is not larger 
than ~ 200 km s _1 . This result is based on the models of Krolik & Kriss (1995) in which 
resonance line scattering is an important process. In this case the X-ray absorber is narrow 
while the UV absorber is broad, implying that the two cannot originate in the same gas. 
We note however that their models covered the b range from 10 km s _1 to 200 km s" 1 , and 
did not definitively rule out the presence of broad (~ 2000 kms -1 ) blended features which 
could mimic a different continuum shape. More importantly, in similar models presented by 
Netzer (1996), absorption and emission due to resonance line scattering almost identically 
cancel out so that no absorption lines are expected even for turbulances significantly larger 
than the thermal line widths. Thus, the non-detection of absorption lines in ASCA data 
may simply be because there are no lines there. Clearly the theoretical situation is too 
uncertain for the absence of features to provide constraints on the b parameter of the 
absorbing gas. Since the UV absorber has a measured width of ~ 2000 km s _1 , we will 
adopt this value in our ensuing discussion of a combined X-ray/UV absorbing region. 

We have seen that the X-ray absorber is consistent with both the presence of the 
broad UV absorption lines in the past and their subsequent absence, but can we develop a 
consistent physical picture to explain this variation? The evolution of ionization parameter 
as seen in figure 4 is suggestive. The earliest known observations of the broad UV absorption 
lines are from 1978. Strong lines were also observed in 1989 (KEA), while they have 
disappeared since ~ 1993. Thus the large scale variability time scale of the absorber seems 
to be some 3-6 years. Given its outflow velocity, the distance traveled by the absorber 
away from the ionizing continuum source in this time is 5 x 10 15 cm. Clearly, this is a small 
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fraction of the distance of the absorber from the central continuum and there would be 
no significant change in ionization parameter corresponding to this change in flux. On the 
other hand, the velocity dispersion of the absorber is large (FWHM~ 2000 km s" 1 ). If the 
absorber is expanding at this rate, its thickness will change by ~ 2 — 4 x 10 16 cm in 3-6 
years. Thus, to have a present thickness of ~ 3.5 x 10 16 cm, the thickness of the absorber 3 
years ago must have been $1.5 x 10 16 cm with a corresponding density of >5 x 10 5 cm -3 . 
Then the ionization parameter in the past could be as low as U=5, (see fig. 4) where 
the strength of the broad UV absorption lines would be in the observed IUE range. CIV 
and NV lines would then be strong. In this scenario the disappearance of the broad UV 
absorption lines is due to increase in U caused by the drop in the density of the absorber 
as it expands while outflowing from the ionizing source. For this scenario to work, the 
fractional change in the thickness of the absorber in 3-6 years should be significant. If the 
absorber is dispersed in velocity space, the 'b' parameter of each cloud would be smaller. If 
we consider a conservative estimate of ~ 200 km s _1 of the velocity dispersion, a significant 
change in the thickness of the absorber would require a density in the higher end of the 
allowed range (n>10 6 cm -3 ). 

The XUV absorber picture is thus completely consistent with the observations: the 
presence of X-ray absorption edges, the presence of broad absorption lines in the old IUE 
spectra and their disappearance in the the new UV observations. In this scenario the broad 
UV absorption lines will not reappear unless a new, similar absorption system is generated. 
Over the next ~ 15 years the current absorber will become more highly ionized, the OVI 
absorption lines will disappear, and the OVIII edge will become stronger than the OVII 
edge. Eventually, it will be completely transparent in X-rays and the OVII/OVIII edges 
will disappear. This is an easily testable prediction with future missions. These estimates 
of time scale and the amount of change in U are based on assuming b=2000 km s -1 and 
density of ~ 2 x 10 5 cm" 3 . If the internal velocity dispersion is smaller, the evolution 
described above would be slower. We note that this picture concerns long term variations, 
we expect in addition small scale variations due to flux variability. 

As noted earlier, new simultaneous HUT (UV) and ASCA (X-ray) observations are 
presented by Kriss et al. (1996a,b). These data show an X-ray ionized absorber similar 
to that reported here and narrow, associated absorption lines in the UV. The authors 
are unable to reconcile the X-ray and UV absorbers with a single absorbing region. This 
is entirely consistent with our model predictions in which the X-ray warm absorber is 
instead associated with the now invisible broad UV absorption. The absorption systems 
in NGC3516 are clearly complex and multiple. The Kriss et al. (1996a,b) study shows 
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that each of the X-ray and narrow UV absorbers require at least two separate systems 
to match the data. More, simultaneous X-ray and UV data allowing monitoring of the 
various components is necessary to disentangle them. We cannot rule out the possibility 
that the X-ray and UV absorbers are different. In that case the disappearance of the broad 
absorption lines may be due to the clouds moving out of line of sight as suggested by Walter 
et al. (1990). Our XUV absorption scenario offers an alternative explanation. 

An edge-on orientation of 3C351, 3C212 and NGC5548 was inferred from their lobe 
dominated radio structure. NGC3516 also has a steep radio spectrum (Ulvastad and Wilson 
1989) typical of lobe dominated AGN. Its radio maps reveal an elongated, one sided, curved 
structure (Miyaji et al. 1992). The geometrical structure of NGC3516 was inferred by Goad 
and Gallagher (1987) by analyzing the velocity field of the circumnuclear emission-line 
region. They found a bipolar outflow from the nucleus in the plane of the sky so that our 
line of sight grazes the putative torus. Thus to date all active galactic nuclei shown to 
contain XUV absorbers are edge on, consistent with the picture of outflowing absorbers 
suggested by Glen, Schmidt, & Foltz (1994). 



7. Conclusions 

ROSAT observations of NGC3516 in 1992 detected an X-ray warm absorber with Nh 
~ 7 x 10 21 cm -2 and U ~ 7.9-12.6. Combining X-ray and UV observations imply that it is 
outflowing with a velocity of ~ 500 km s" 1 and has an internal velocity dispersion of <2000 
km s -1 . The constraints on density are not very tight: 10 3 <n<10 7 cm" 3 . The distance 
of the absorber from the central continuum source is <5 x 10 18 cm and its thickness is 
about 3.5 x 10 15 cm assuming density of 2 x 10 5 cm" 3 . The physical characteristics of the 
absorber are consistent with all the observations: the presence of X-ray absorption edges, 
the presence of broad absorption lines in the old IUE spectra and their disappearance in 
the the new UV observations. In this scenario of the XUV absorber in NGC3516, the 
disappearance of the broad UV absorption lines is due to the current high value of U caused 
by the drop in the density of the absorber as it expands and moves away from the ionizing 
source. If this dynamical model is correct, we expect the OVII absorption edge in the X-rays 
to weaken compared to the OVIII edge. Eventually the source would be transparent to the 
X-rays as the absorber continues to expand, unless a new absorption system is produced. 
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Figure 1: NGC 3561 light curve (top panel), t vii (middle panel), t viii (bottom 
panel), for the eight GTIs. 

Figure 2: Residuals to the simple power-law fit to the entire data. The strong negative 
residuals between 0.7 and 1 keV are clearly visible. Channels 13-25 corresponding to this 
energy range were ignored in the PL fit to show the absorption edges clearly. 

Figure 3: Data and the best fit spectrum (PL+2 Edges, N# fixed to Galactic). 

Figure 4: Ionization fractions f of OVI, OVII, OVIII, CIV and NV as a function of 
ionization parameter, U. The vertical lines define the range of U for which the ratio 
iovn/^ovin hes within the observed range. The arrows on the CIV and NV curves indicate 
the lower limits of f C /y^3 x 10" 4 and i NV ~3.1 x 10" 4 based on the published IUE data 
(see text). 

Figure 5: The input spectrum (solid line) and the transmitted spectrum (dotted line) for 
N# = 10 22 cm" 2 and U=10. The transmitted spectrum shows strong OVII/OVIII edges, 
but no Fe-K edge or low energy cutoff. 
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Table 1: ROSAT Observations of NGC 3516 



Observation 


Net Counts 


Exposure 

(s) 


Net Count Rate 


Total 


59225±246 


13081 


4.53±0.02 


GTI 1 


12860±127 


3176 


4.05±0.04 


GTI 2 


4186±69 


1150 


3.64±0.06 


GTI 3 


4068±71 


1014 


4.01±0.07 


GTI 4 


14848±116 


2900 


5.12±0.04 


GTI 5 


7007±96 


1375 


5.10±0.07 


GTI 6 


7670±97 


1612 


4.76±0.06 


GTI 7 


1027±39 


232 


4.43±0.17 


GTI 8 


7571±97 


1622 


4.67±0.06 
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Table 2: Spectral fits to ROSAT data of NGC 3516 



Data 


Model 


a E 


Normalization 0- 


T OVII 


TOVIII 


X 2 (dof) 


F b 


GTI r 


PL 

PL+ lEdge 
PL+ 2Edges 


1.22±0.01 
1.07±0.03 
1.04±0.035 


1.36±0.01 
1.70±0.07 
1.76±0.07 


0.8±0.15 
0.46«i? 


0.45±0.19 


132.2 (30) 
43.55 (29) 
31.11 (28) 


59.0 
11.2 


GTI 2 


PL 

PL+ lEdge 
PL+ 2Edges 


1.18±0.04 
1.06±0.06 
1.05±0.06 


1.24±0.04 
1.51±0.11 
1.53±0.12 


o.7itg£ 

0.62±0.35 


n 1 1 +0 - 36 


56.5 (30) 
33.68 (29) 
33.45 (28) 


19.6 
0.2 


GTI 3 C 


PL 

PL+ lEdge 
PL+ 2Edges 


1.23±0.03 
1.08±0.06 
1.02±0.07 


1.33±0.03 
1.69±0.13 
1.8±0.14 


0.9±0.3 
0.34±0.33 


u - ' °-0.37 


62.87 (30) 
32.17 (29) 
22.00 (28) 


27.7 
12.9 


GTI 4 C 


PL 

PL+ lEdge 
PL+ 2Edges 
Best Fit d 


1.26±0.01 
1.12±0.02 
1.09±0.03 
1.24±0.09 


1.69±0.02 
2.08±0.05 
2.16±0.09 
2.23±0.09 


0.77±0.08 

o.43l8:ii 

u - '°-0.25 


0.45±0.18 

o.i9±8:?i 


154.1 (30) 
59.57 (29) 
44.87 (28) 
34.55 (27) 


46.0 
9.2 
8.1 


GTI 5 


PL 

PL+ lEdge 
PL+ 2Edges 


1.22±0.02 
1.09±0.04 
1.08±0.05 


1.72±0.08 
2.09±0.11 
2.13±0.13 


n 79+0.20 

u - ' z -0.19 

0.571^ 


0.201°;^ 


66.20 (30) 
25.16 (29) 
23.91 (28) 


47.3 
1.5 


GTI 6 


PL 

PL+ lEdge 
PL+ 2Edges 


1.16±0.02 
1.07±0.04 
1.06±0.04 


1.55±0.03 
2.1±0.1 
2.1±0.1 


1.1±0.2 
1.0±0.3 


0.06jl8;g 6 


117.4 (30) 
20.89 (29) 
20.81 (28) 


134.0 
0.1 
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Table 2: Continued. 



GTI 7 


PL 

PL+ lEdge 
PL+ 2Edges 


1.2±0.1 
1.0±0.13 
l.OiO.l 


1.5±0.1 
2.0±0.3 
1.99±0.35 


1 o +0 - 7 

n 7Q+0.36 
u - ' °-0.68 


o.3l8:I 


19.91 (30) 
11.83 (29) 
11.5 (28) 


19.8 
0.8 


GTI 8 C 


PL 

PL+ lEdge 
PL+ 2Edges 


1.19±0.02 
1.06±0.04 
1.03±.0.05 


1.60±0.03 

1.9±0.1 
2.01±0.11 


0.7±0.2 

n 4+ - 25 


0.4+^ 


71.16 (30) 
29.08 (29) 
22.74 (28) 


42.0 
7.8 


GTI 1-8 


PL+ 2Edges 
PL+ 2Edges 
(N H free) 


1.06±0.02 
1.16±0.04 




0.53±0.09 

o.74±8:il 


0.37±0.09 


244.6 (245) 
226.9 (244) 





a: in units of 10~ 2 photons keV -1 cm~ 2 s _1 at 1 keV 

b: Parameter of the F-test. Note that F > 7.5 is 99% significant and F > 13 is 99.9% significant 
c: datasets for which two edges are required by the data 
d: PL+ 2Edges+ Excess absorption at the source (see text). 
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Table 3: Warm absorber model fits to ROSAT data of NGC 3516 



Data 


OLE 


Normalization" 


N# 6 (warm) 


log 17 


X 2 (dof) 


GTI 1 


o.9il8:8i 
1 1+°- 2 

1-1-0.1 


9 1 cr+0.23 
z - io -0.17 
9 99+O.I5 
z - zz -0.17 


i.26l8:^ 
o.77l8:| 


i.34l8:S 
i.09l8:I 5 


37.02 
30.41 


(28) 
(27) 


GTI 2 


0.95±0.11 
1.4±0.3 


1 7O+0-28 
l-<y_0.22 
9 9+O.6 
z - z -0.4 


0.69±0.37 

0.418:? 


i.i6l8:l 
o.45l8:^ 


32.56 
26.41 


(28) 
(27) 


GTI 3 


0.82±0.65 

0.8218:2 


9 7+0.48 
z -°'-0.03 

2.41181 


3.19±0.85 

q qpr+0.79 
°-° d -3.34 


1.4918:8^ 
i.49l8:88 8 


23.80 
23.52 


(28) 
(27) 


GTI 4 


o.9±8:f 
1.318:? 


2.76±0.12 

2.918:1 


2 1+ L2 

z - i -2.05 

0.618:1 


1.4618:8! 
o.9l8:i 


50.72 
31.54 


(28) 
(27) 


GTI 5 


0.94±0.56 

i q+0.5 


2-593:18 

2.81I8I 


1.2111 

0.5I8I 


1-341818 

0.81I8I 


26.94 
17.35 


(28) 
(27) 


GTT fi 


W . U xj _1_ W . t: (J 

i.3l8: 2 7 


2.9l J 


90+9,-L 

o.6l8:i 


1 14+0-19 
1-1^-0.18 

0.8I8I 


32.00 
21.45 


(28) 
(27) 


GTI 7 


o.87l8:| 

l-Olo^ 


2.44±8:H 

9 c+2.5 
z -°-0.7 


1.0315-g 4 

o.8l8: 9 7 


1 99+O.27 

1.06l?1 6 


11.98 
11.77 


(28) 
(27) 


GTI 8 


0.86±0.49 
1 1+0.2 

i - i -0.3 


2.57±0.43 

9 fi +0.2 
z -°-0.3 


i.93li:S 


1-4518:? 4 
1 1+0.3 

1.1—0.4 


20.71 
19.05 


(28) 
(27) 


Total 


1.21±0.06 


z - J ^-0.08 


0.7±0.1 


1.0±0.1 


28.01 (27) 



a: in units of 10 2 photons keV 1 cm 2 s 1 at 1 keV 



b: in units of 10 22 atoms cm 2 
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Fig. 3.- 
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